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ABSTRACT 


Nonstationary vibratory phenomenon is exhibited by machines and components, when 
they are subjected to variable frequency operations like start-up and shut-down. While 
nonstationary behaviour of linear systems are fairly well documented, information on 
nonlinear systems is scant due to lack of analytical procedures. The present study 
attempts to carry out numerical and experimental exercises to understand such behaviour. 
The exercise is restricted to systems involving spring nonlinearities only. Damping is 
taken to be linear in the entire analysis. Four different types of spring nolineeirities are 
studied - (i) polynomial nonlinearity; (ii) symmetric square nonlinearity (iii) bearing 
nonlinearity (iv) bilinear nonlinearity. Single degree freedom systems are considered and 
their response is numerically obtained through fourth order Rxmge-Kutta technique, 
employing MATLAB Simulink Toolbox. The results in the first case are compared to 
those obtained using the approximate asymptotic analysis. In the remaining cases system 
response is sought to be understood through equivalent piecewise polynomial form 
expressions. Experimental investigations are carried out through Virtual Instrumentation 
techniques for signal generation and data acquisition. The response of a cantilever beam, 
subjected to frequency swept excitation is reported. 



CHAPTER 1 


INTRODUCTION 

In most vibration analysis problems an acceptable idealization and simplification is made 
- that the system components are time independent and linear. Thus masses, restoring 
forces, frequencies and amplitudes of external excitation are considered to be time 
independent. Oscillatory motion of such a system is called stationary. When one or more 
of the system or external excitation components vary with time, the system is said to be 
non-stationary and the oscillatory motion of such a system is non-stationary. 
Nonstationary vibratory phenomenon is exhibited by machines and components, when 
they Eire subjected to variable frequency operations. Start-up and shut-down operations of 
machinery are common instances of occurrences of such phenomenon. Many 
reciprocating engines, turbines are operated above their critical speed and pass through 
this critical speed during start-up and shut-down processes. The dynamic stresses 
experienced during such processes and their duration need to be analysed for safe 
operations. The objective of the present study is to analyse nonstationary vibrations of 
nonlinear systems under variable frequency excitation. While procedures for analysis of 
linear systems subjected to excitation forces, with frequency varying with time, are fairly 
well established, very few studies are available on the behaviour of nonlinear systems 
under such excitation conditions. For linear systems, it is known that the maximum 
amplitude during a frequency sweep will be smaller than if the excitation J&equency was 
held constant, and further, that an apparent shift in the position of critical will occur. The 
lack of information of such nonstationary phenomenon of nonlinear systems, is primiirily 
due to the considerable mathematical complexity of such problems and a complete 
absence of closed form analytical solutions. 


The problem of non-stationary oscillations in mechanical and electrical systems was 
discussed by Lewis (1932). He presented an exact solution of the problem of running a 
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system having a single degree of freedom and linear damping through its critical speed 
from rest at a uniform acceleration rate. The analytical solutions of non-stationary 
vibrations were first proposed by Krylov and Bogolyubov (1961), through their method 
of asymptotic expansions. They gave an effective method of investigation, i^plicable to 
an extensive class of oscillatory processes described by nonlinear equations with varying 
coefficients. Reference can be made to studies by Gasch et al (1979), Vyas et al (1986), 
Irretier (1984), Nonami et al (1988) for analysis of problems concerning rotating shaft 
and blading systems during acceleration through critical speeds. While the stationary 
response of nonlinear systems is characterised (as distinct from stationary response of 
linear systems) by presence of subharmonic frequencies and an inclined backbone curve, 
information on the response during nonstationary conditions is scant. An asymptotic 
approximation of the envelope of the response of a Duffing oscillator under linearly 
varying excitation frequency, based on slowly varying parameters was developed by 
Mitropol’skii (1965). Agrawal and Ewan-Iwanowski (1976) investigated non-stationary 
response for various kinds of summed-and-differential harmonic oscillations. They 
obtained approximate solutions by means of the asymptotic method and gave amplitude 
variation curves calculated numerically. In experiments were also carried out for 
verification of analytical results. Recently, Ishida et al (1992) studied non-stationary 
oscillations of rotating shaft with nonlinear spring characteristics during acceleration 
through a major critical speed. Approximate solutions for steady state and non-stationary 
oscillations by asymptotic method and proposed a complex FFT method. The method is 
based on the signal processing concept. Spectral distribution is obtained for the entire 
time history of accelerated or decelerated signal. A desired fi^uency component is 
extracted by filtering. This extracted frequency component is transformed into time 
history by inverse FFT algorithm. The amplitude variation of this transformed signal is 
then plotted. 

A major limitation of the studies on nonlinear systems is that procedures of analysis exist 
only for a very limited class of nonlinear functions. Moreover the solutions are 
approximate and apply with reasonable accuracy only within certain range of operations. 
If the amount of non-linearity is not too large, analytical methods may be used to yield an 
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approximate solution. If the amount of non-linearity is larger, analytical methods may not 
be sufllcient and solution may be possible only by numerical or graphical methods. Thus, 
any solution applies for only one particular set of conditions. Furthermore, the process of 
obtaining solution is usually tedious, requires considerable manipulation if good accuracy 
is to be obtained. The focus of the present study is to study the nonstationary response 
behavior of linear and nonlinear single degree of freedom systems, through numerical 
simulation. Tools available in MATLAB software have been employed for numerical 
simulation. Various kinds of commonly occurring nonlinearities in mechanical systems 
have been investigated. The analysis is carried out in terms of non-dimensional 
parameters, to the extent possible; The response of linear systems subjected to variable 
frequency excitation is presented in Chapter 2. Response of four different types of 
nonlinear systems is investigated in Chapter 3. The results are analysed as ftmetions of 
nonlinear parameters, damping and excitation frequency acceleration rates. Experimental 
work carried out on a laboratory set-up is described in Chapter 4. Conclusions and scope 
for further work are given in Chapter 5. 
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CHAPTER 2 


NON-STATIONARY VIBRATIONS OF LINEAR SYSTEMS 

Response of a linear single-degree of freedom system under frequency swept excitation is 
presented in this Chapter. This analysis is readily available in literature and is presented 
here, for the sake of completeness and later comparison with the response of nonlinear 
systems. 

2.1 Equation of motion 

For a simple mass, spring, dashpot system. Figure (2.1), acted upon by a sinusoidal force, 
with initial frequency (Dq and changing at a constant rate a , the equation of motion can 
be written as 


mx + cx + kx = fQSin{cjt)J + -at^) 

2 


( 2 . 1 ) 


where a dot denotes differentiation with respect to time t, k and c are constants of 
proportionality for restoring force and damping respectively, m is the mass of the system. 


Defining the following non-dimensional parameters 

T = pt, r^^oyjp, r^=alp\ 

x=xfx^i, 2^ = clk p = ^klm 

with =fjk 


equation (2.1) is rendered in the following non-dimensional form: 

x" + 2^' + x = s\R{r^r + ^r^v^) (2.2) 
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2.2 Response Simulation 

The response of the nondimensional governing equation (2.2) is simulated using 
MATLAB (Simulink toolbox). A sample program for the simulation of equation (2.2) on 
MATLAB is shown in Appendix 1 . The numerical simulation is based on fourth-order 
Runge-Kutta technique. 

Fig. 2.2 shows a typical numerically simulated response of equation (2.2), for a damping 
ratio, ^ = 0.01 and rate of acceleration = 0.001 . The sweep is carried out from initial 

frequency, = 0 to a final nondimensional frequency value equal to 2.5. The response is 

shown by solid line and the dark solid line shows the envelope of the response, which 
defines the maximum amplitude attainable, during the process. Henceforth, in the present 
study, response is displayed only through the envelope curves, for better understanding of 
qualitative behavior of the system. 

The response is studied as a function of different rates of acceleration and damping ratios. 
Fig. 2.3 describes the response envelope for three frequency sweep rates = 0.001, 
0.003 and 0.005. The damping ratio is kept constant (^ = 0.01 ). The stationary frequency 
response curve = 0 , is also shovm in the figure. The resonant response amplitude, can 

be readily seen from the stationary response curve, to be equal to 50, which as expected is 
equal to the quality factor 1 / 2^ . The resonant response frequency, on the 

nondimensional instantaneous frequency {r = az) axis can be seen to be equal to 1, 
For the nonstationary cases ( ^ 0 ), the response can be seen to attain an apparent 

resonance, with maximum amplitude less than the resonant amplitude of the statipnary 
case. Also this apparent resonance occurs not at instantaneous excitation frequency value 
r= 1, but at an instant shifted to the right of r = 1. These apparent resonant amplitudes 
can be seen to decrease with increasing rates of acceleration. Also the shift away from the 
instant r = 1 , becomes more pronounced with increase in the acceleration rate. This can 
be explained by the fact that while the phenomenon of resonance does occur when the 
excitation frequency is equal to the system’s natural frequency (i.e. r =1), resonant 
amplitudes also need time to build up. 
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Non-Dimensional displacement ^ 






Fig 2.5 Acceleration and Deceleration for linear system 


In the case of stationary vibrations at a constant excitation frequency, the system reaches 
its steady state value. However, sufficient time is not available to the system to reach 
steady state, when it is subjected to a frequency sweep. The amplitude of vibration does 
start increasing, when the instantaneous excitation frequency momentarily becomes equal 
to the natural frequency. In addition to the externally provided energy being utilised in 
building up resonant amplitudes, a portion of this energy is also utilised in accelerating 
the system. Under the influence of these two activities an apparent resonance is exhibited, 
with a vibratory amplitude, less than what would have been achieved, if the process were 
stationary, and at an instant shifted in the direction of frequency sweep. The higher is the 
acceleration rate, less is the time available to build resonant amplitudes, as seen in Figure 
2.3. In absence, of a steady state, system can also be seen to exhibit the beat 
phenomenon, in the vicinity of the region r = 1 . Fig. 2.4 shows the influence of the 
damping const , and as expected higher values of damping effect greater suppression of 
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the resonant amplitude, without any other qualitative change in the behaviour. The 
behaviour under a negative frequency sweep rate, is shown along with an identical 
positive sweep rate in Fig. 2.5. The response pattern is practieally symmetric about r = 1, 
with the maximum amplitude being attained at an instant shilled to the left of r = 1, 
which in the case of negative sweep rate is the direction of frequency change. The initial 
frequency in the case of negative sweep is 2.5, while the curve for positive sweep rate 
starts with a value of r-g = 0 • The minor difference in the peak response values is due to 
the difference in the initial conditions, in the two cases. 
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CHAPTER 3 


NONSTATIONARY OSCILLATIONS OF 
NONLINEAR SYSTEMS 

Nonlinearities in mechanical systems can arise from several sources. These are broadly 
classified into (i) geometrical nonlinearities - arising due to relatively large 
displacements and (ii) physical nonlinearities - arising due to nonlinear or multi-valued, 
stress-strain relationships. This non-linearity may be present due to stiffening 
characteristics of material itself. For instance, an organic material such as rubber or 
leather used in couplings and vibration absorbers, have their modulus of elasticity 
increasing with elongation. Another example of variable flexibility is encountered in the 
case of cast iron or concrete structures. Elastic characteristics can be seen to vary with 
displacements in some steel springs as well. Natural frequencies of systems involving 
such elements depend on the amplitude of vibration. Another kind of nonlinearity is 
encountered in systems where the damping forces caiuiot be simply represented by linear 
functions of velocity. For instance, resistance of air or liquid, at considerable high speed, 
is taken to be proportional to the square of the velocity. The present study is restricted to 
systems involving spring nonlinearities only. Damping is taken to be linear in the entire 
analysis. Four different types of spring nolinearities are studied - (i) polynomial 
nonlinearity; (ii) symmetric square nonlinearity (iii) bearing nonlinearity (iv) bilinear 
nonlinearity. Single degree freedom systems are considered and their response is 
numerically obtained through fourth order Runge-Kutta technique, employing MATLAB 
Simulink Toolbox. The results in the first case are compared to those obtained using the 
approximate asymptotic analysis for the envelopes, as developed by Mitropol’skii (1965). 
In the remaining cases equivalent polynomial form expressions are obtained for the 
nonlinear functions for the interpretation of the numerical results. Results are obtained for 
both positive and negative sweep rates. The analysis is carried out in nondimensional 
form and the significant features of the nonstationary response associated with 
nonlinearity forms under consideration are presented. 


10 



3.1 Polynomial nonlinearity 

A single degree of freedom system, with linear damping, excited by a sinusoidal force 
with a frequency constantly increasing with time, as given below is considered. 


mx + cx + f{x) = /o sin(iyQi 


(3.1) 


fix) is the non-linear restoring force, which is taken to be of the following general 
polynomial form. 


fix) = k^x -f + 


(3.2) 


Restricting to terms up to cubic, the following two cases are investigated 
Case (1) Cubic nonlinearity alone i.e. , k 2 =Q, . 

Case (2) Combined Square and Cubic nonlinearity i.e. A:, 0 , k2^0, kj^O. 


3.1.1 Case (11 Cubic nonlinearity alone ( A] 0 . = Q • ^3 ^ Q ^ 

A system with cubic nonlinearity is the commonly known Duffing oscillator and is 
considered to be a benchmark system in nonlinear analysis. The governing equation (3.1) 
becomes 


mx-i-cx-i- A:,x + AijX^ = /q sin(mot-i--^a/^) 


(3.3) 


Response is numerically computed for both hard ()9>0) and soft(/0<O) springs. 
Defining the following parameters 

^ = ro=cojp, r„=a/p^ 

x = x/x,„, x' = x'/x,,, x'’ = xfx^ (3.3a) 

^ = 2C = clmp p^^jkjm, 

withx„ =/o/A:, 

equation (3.3) can be written in non-dimensional form as, 
x" -I- 2^' + X + .^x^ =sin(roT + '^r„T^) 


(3.4) 



In the above ^ is the linear damping ratio and A is a nondimensional nonlinear 
parEimeter. It is to be noted here that A , in addition to the nonlinear stiffness term ^ 3 , 
involves the excitation force amplitude /q also. The response of equation (3.4) is 

numerically simulated in MATLAB. The Simulink model developed for the purpose is 
given in Appendix 2. The nondimensional response x is plotted in Fig. 3.1 (a) as a 
function of the instantaneous nondimensional excitation frequency r(=ro+r^T), for 
nonlinear parameter A = 0.05 . The plot for the stationary response of the system is also 
shown dashed in this figure. This plot is obtained from the steady state response of the 
system for constant frequency excitation, as given by the following governing equation 
3c'' + 243c' + 3c + A3c^ = sin(rr) (3.5) 


( r is a constant, in the above) 

The stationary response curve shows the features like the backbone curve and the jump 
phenomenon, characteristic to the Duffing oscillator. In contrast to a linear system the 
natural frequency of a Duffing oscillator is amplitude dependent and is characterised by 
the backbone curve. The backbone curve describing the relationship between the 
nondimensional nonlinear parameter A , nondimensional excitation frequency r and the 
nondimensional amplitude A of vibration, can be written as (Stoker, 1950). 


/l(l-r') + 3A— +( 2 ^r)^=l 

4 


A two term Volterra series approximation of the response of equation (3.5), for small A 
can be vmtten as (Schetzen, 1979) 

x(r)«Re{//,(;>>^"} + ^Re{/73(;>,7>.-yr)e^^^} + ^Re{//3(7>,yr,yr)e^^"} 

4 4 

= Re{//, Ur)e ^'^ } + 1 A Re{//,' (y>)i/,* (y>)e^" } + jRe{Aif/ (yr)//, (3yr)a^^" 


x(T) = \H^ (;>)|cos[rr-(9(r)]+3-|/f,(yr)|'' cos[rr-2^(r)] 
+ f l«i 0>)r |W, 0-3r)| cos[3rr - 30(r) - S(3r)] 
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Non-Dimensional instantaneous frequency r(= +r^r) 

Fig. 3.1(a) Steady and non-stationary response of cubic non-linear 

hard system (A > 0 ) 



Non-Dimensional instantaneous frequency r(= -i-r^r) 

Fig. 3.1(b) Steady and non-stationary response of cubic non-linear 

soft 8ystem(/l. < 0) 


B 




Non-Dimensional displacement x Non-Dimensional displacement x 



0 0.5 1 1.5 2 2.5 3 3.5 4 4.5 5 


Non-Dimensional instantaneous iBrequency r(= +r^r) 



Non-Dimensional instantaneous frequency r(= + r^r) 


Fig. 3.2(a) Influence of cubic non-linear parameter A during acceleration 

(hard spring, A >0) 
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Non-Dimensional displacement x 



Non-Dimensional instantaneous frequency r(= + r^r) 


Fig 3.2(b) Influence of cubic non-linear parameter A during deceleration 

(hard spring, A >0) 
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In the above expressions //, {jr) is the first order Frequency Response Function (FRF) 
of the equation and is given by 

H,Ur) = \m-r^) + mr)\ (3.7a) 

The response given by equation (3.7) comprises of harmonics at the excitation 
frequencies r and a subharmonic at firequency equal to 3r . The amplitude of the first 
harmonic then can be written as 

A =|H,(;r)l 1+3^1//, (;>)|%-^^('^ (3.8) • 

and that of the subharmonic as 

(3-9) 

The subharmonic resonance can be seen in Fig.3.1(a). 

The nonstationary response in Fig.3.1 (a) has been computed for both positive and 
negative rates of acceleration, . The magnitude of is 0.001. During positive 

acceleration, the initial frequency is taken as zero, while in the negative acceleration 
case the initial frequency, is equal to 5.0 and the frequency range is swept back to 

zero. The corresponding behaviour for a soft spring (A = - 0.01) is shown in Fig.3.1 (b). 
Fig. 3.2 (a) shows the influence of various values of the nonlinear parameter A , for hard 
spring case (A positive) with acceleration rate =0.001, while Fig 3.2(b) shows the 

corresponding behaviour for a negative acceleration. Results for various negative values 
of A are shown in Figs.3.3 (a), (b). The influence of different acceleration and 
deceleration rates for hard spring is depicted in Fig. 3.4. The influence of damping on the 
overall response is shown in Fig.3.5 with acceleration rate = 0.001. In Figs. 3.2 to 3.5, 
wherever possible, subharmonics are zoomed and presented. 
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Non-Dimensional displacement x 



Fig. 3.3(a) Influence of cubic non-linear parameter JL during acceleration 

(soft spring, A<0) 







Non-Dimensional instantaneous jfrequency r(= + r„r) 



Non-Dimensional instantaneous frequency r(= + r^r) 


Fig. 3.3(b) Influence of cubic non-linear parameter X during deceleration 

(soft spring, X <0) 
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Lcement x 


A = 0.05 


Non-Dimensional instantaneous frequency r(= + r^r) 
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Non-Dimensional instantaneous frequency r(= r,, + r^r) 

\ Influence of acceleration rates (hard spring, A > 0 ) 













The following conclusions can be drawn from Figs. 3.1 to 3.5 on the nonstationary 
vibrations of the Duffing oscillator - 

(a) The maximum response achieved by the system decreases with increiise in the 
acceleration rate. The sharpness of the resonance curve also becomes less pronounced 
as the transition through resonance becomes faster. This behaviour is similar to that of 
a linear system. 

(b) The resonance curves are not symmetric as in linear systems. The response during 
positive acceleration is starkly different than that during negative acceleration. 

(c) The nonstationary response of a nonlinear system is associated with a jump 
phenomenon. The jump is steeper, when the rate of acceleration is less. 

(d) Higher nonlinearity results in lower values of the maximum response achieved during 
positive acceleration The shift in the apparent resonant point is higher for higher 
magnitudes of the nonlinearity parameter A , during positive acceleration, while for 
negative acceleration the apparent shift of the resonant point towards left is less for 
higher values of A . This is in accordance with the fact that for higher nonlinearity, 
the bent of the backbone curve to the right of r = 1, is more pronounced. 

(e) A phenomenon exactly opposite to that described in (d) above, takes place when the 
spring is soft (negative A ). 

(f) The nonstationary response is also characterised by the presence of transient 1/3 
subharmonic. The subharmonic also undergoes an apparent shift, depending on the 
magnitude and sign of nonlinearity as well as the magnitude and sign of the 
acceleration rate. 

(g) The influence of linear damping, as expected, is a reduction in vibratory amplitude. 
Comparison with Asymptotic Method 

The numerically simulated results have been compared with those obtained through the 
Asymptotic Method proposed by Krylov and Bogolyubov (1956) through their method of 
asymptotic expansions. This has been done in order to check the correctness of the 
numerical simulations. A brief description of the asymptotic method is given below. 
Writing a general governing equation for a non-linear single degree of freedom system as 

x-¥ p^x = (3.10) 
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where p = yjk/m is the linear natural frequency and £• is a small parameter, the purely 
harmonic response for the case when £■ = 0, is 

x = acosy/ (3.11) 

with a constant amplitude a and uniform rotating phase angles , i.e. 

a = Q,\i/ = <a (3.12) 

The existence of nonlinear perturbation ( £ ?* 0 ), results in the appearance of overtones in 
the solution of equation (3.10). It also establishes dependence between the instantaneous 
frequency y/ and the amplitude a, and finally gives rise to systematic increase or 
decrease in the amplitude of the vibrations, depending upon whether the energy is 
expelled or absorbed by the perturbing force. Taking all these into account a general 
solution is sought in the form of expansion 

x = acos(ii/)+eu^{a,\fr)^£'^u^{a,\f/) + ....-\-£'”u^{a,y/) m = l,2,... (3.13) 

where u^{a,\i/),U 2 {a,\i/) are periodic functions of the angle ij/ with period 2k and the 
quantities as function of time, are defined by the differential equations 
a = £A^{a)+ £^ A-^ia) +..,*£'” A„{a) (3.14) 

= Q) + £B^{a)+ £^ B2{a) +... + £"' B^{a) 

Thus, the problem of constructing asymptotic approximations for equation (3.10) reduces 
to the choice of suitable expressions for the functions 
u^(a,y/),U 2 (a,y/) ...,A^{a),A 2 ia)...,B^{a),B 2 ia)... in such a way that equation (3.13) 
would, after replacing a and \i/ by the functions of time defined by equation (3.14), 
serve as a solution of the initial equation (3.10). The practical applicability of asymptotic 
method is known to be determined not by the convergence properties of the series (3.13) 
and (3.14) when m^oo, but by their asymptotic properties for a given fixed value of m 
when £ -> 0 . Applying the above method to the cubic non-linear case under present 
consideration (equation reproduced below) 


22 



mx + cx + k^x + sin((y5r +-^crf^) 

and further rendering it in nondimensional form as 

ic" + 243c' + jc + 3 c^ =;'sin(roT + ^r„r^) 
with 

T = pt 

x = {jkjr^)x 

2^ = cl mp 


ro=o)Jp 

= {^kjk, y 


3F' 


r = j^,«AAi • 

the first approximation of the response is 

X = acos(0 + ^) 
with 


fa = Clip'' 

^' = (7*3 At V 
p = ylkjm 


a = -^a - 


r 


\ + v 


a 


■cos 


(«>) 


<f> = l-v + 3 — + 


■5\n{<f>) 


8 a(l + v) 

The second approximation can be shown to be (Evan-Iwanowski ,1976) 

X = acos(^ + ^)+-^cos3(0 + ^) 


with 


• ^ 3aV ( 

a = -La H — - Y 

8 


2a^{l-v) 

l + y 8(3-vX1 + k)^ 


cos(^) 


+r 


1 


V-- 




(l+vX (1 + v) 


sin(^) 


j , 3a' ISa' , r 

^ 8 2 256 a 


1 3a'(5-3v) 1 


/ 

+ 

V V 


V- 


c 


(i+v)' (i+v y 


cos(^) 


(3.15) 


(3.16) 


(3.17) 

(3.18) 


(3.19) 


(3.20) 
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Fig. 3.6 Response comparison of simulation and asymptotic method 


The results with the second approximation (equations 3.19, 3.20) for 

= 0.001 = 0.01 , and y = 0.03 are shown in Fig 3.6. The numerical results obtained 

through the fourth order Runge-Kutta model in MATLAB are also shown in the figure. 
The good qualitative and quantitative closeness between the two sets of curves is taken as 
a measure of the correctness of the numerical model developed in MATLAB. 

3.1.2 Case r21 : Combined Squa re and Cubic nonlinearitvt 

A system with both square and cubic stiffiiess nonlinearity is considered next The 
governing equation with linear damping and a sinusoidal excitation force with firequency 
constantly increasing with time can be written as, 

mx + cx + k^x + k.^x^ +kjX^ = sin(<j)o/ + ^at^) (3.21) 
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Fig. 3.7(a) Non-stationary response of square-cubic and cubic nonlinear systems 



Non-Dimensional instantaneous frequency r(= +r„r) 


Fig. 3.7(b) Subharmonics for combined square-cubic and cubic nonlinearity 
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Non-Dimensional instantaneous frequency r(= Tq + r^x) 

Fig. 3.7(c) Non-stationary response of square-cubic (/9<0,A>0) 
and cubic nonlinear systems 


Employing the parameters defined in equation (3.13 a) and ^ = ^ 2 / 0 /^!^ » equation 
(3.21) can be written in nondimensional form as, 

x” + 2^' + x-^ +2jc^ =sin(r(,r + -^r„r^) (3.22) 

The MATLAB Simulink model for numeric solution of equation (3.22) is given in 
Appendix 3. The nondimensional response of the system is shown in Fig 3.7(a), for a 
nonstationary case, with acceleration, = 0.001 . It can be seen that the presence of 

square nonlinearity, causes r/2 subharmonic resonance. This subharmonic has been 
zoomed in Fig. 3.7 (b). Other characteristics appear similar to a cubic nonlinear system. 
Fig 3.7 (a) corresponds to a positive value of (=0.025), while in Fig 3.7(c) the response 
corresponds to a negative value of (= -0.025). In both cases the value of cubic 
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nonlinear parameter A and rate of acceleration were kept constant at 0.05 and 0.001 

respectively. The figures also show comparison with the response of a purely cubic 
system {(5 =0). It can be noted from Figs. 3.7 a, c that a positive value of , results in 
downward shift of the response curve in comparison with that of a purely cubic system. A 
negative value of /? causes the opposite. ( This fact can be understood by transforming 

equation (3.21) into an alternative form mq + cq-^ a^q + a^q^ = sm{a)Qt + a/ V2) + ; 

with q = x-\- being constants depending upon and ; and carrying 

out further analysis as suggested by Hayashi, (1964). 

3.2 Symmetric Square Nonlinearity 

A system with nonlinear restoring force g(x)=xlxl is considered here. The governing 
equation is written as 

nvc + cx + h:\x\= + ^ca^) (3.23) 

and in nondimensional form as 

x" + 2^' + X I X 1= sin(ror + ^r^r^) (3.24) 

where =kXf^jm\ = f^/k (other parameters same as defined in equation 3.13 a) 
The Simulink model developed for the purpose is shown in Appendix 4. Fig 3.8 shows 
the response characteristics of equation (3.24) for acceleration and deceleration rate r„ = 
0.001. The response is characterised by the presence of both 1/3 and 1/2 subharmonics. 
Explanation of this behaviour can be attempted by expressing the restoring force 
g(x) = X I X I , by an equivalent polynomial power series form of the type given by 

equation (3.2) i.e. g(x) = + + It is to be noted here, that the 

equivalent terms k^,k 2 ,kj will be dependent on the response level x. If the curve fitting is 
carried out over the entire range of the maximum response level (-14 ^ x ^ 14) of Fig.3.8 
the following values of equivalent terms are obtained, 
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Fig. 3.8 Response of symmetric square system during acceleration 

and deceleration 


it, =4.3904 

itj =0.0 (3-25) 

it3=0.0519 

Fig. 3.9 gives a comparison of the actual square symmetric spring force and that obtained 
by the equivalent polynomial representation. Fig. 3.10 shows a comparison between the 
behaviour of the system with actual nonlinear function g{x) = x\x\ and that with 

equivalent polynomial form function g{x') = k^x + k 2 X^ + k^x , with values as given 
above in equation (3.25). The curves have been obtained for acceleration rate 
= 0.001, = 0.01 .While the equivalent polynomial approximates the symmetric square 

nonlinearity closely in the resonance zone, it misses out on depicting the subharmonics. 
However, if the curve fitting to obtain the equivalent polynomial form is restricted to the 
response 
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Fig. 3.10 Comparison of square symmetric jc | x | and equivalent 
polynomial system response 
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Non-Dimensional instantaneous frequency r(= +r^r) 


Fig. 3.11 Subharmonics of square symmetric x\x \ system and equivalent 

polynomial form 
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level in the subhamionic zone (—2 ^ x £ 2) the values of coefficients of equivalent power 
series, come out to be, 

A, =0.6265 
= 0.0 

yfc, =0.3638 

The response behaviour for this polynomial approximation (Fig.3.1 1) can now be seen to 
reflect the subharmonic response characteristics of the x | x | system, quite correctly. 

3.3 Bearing Nonlinearity 

Most rotors are supported in oil-film or in rolling element bearings. The bearings 
influence the rotor vibrations to greater or less degree by their dynamic properties in 
comparison to those of the rotor and the remaining components of the assembly. This 
influence originates essentially from the ratio of rotor stiffiiess to bearing stiffness. For a 
relatively elastic rotor the influence of the bearings is small while dynamics of a stiff 
rotor are determined by the bearings. A typical case of a rigid rotor in rolling element 
bearings with sinusoidal exciting force accelerating linearly with time is considered. The 
force-deflection relationship for a ball-bearing is generally expressed as, (Kramer, 1993) 

g(x) = kx^ 

= ib:|x|^ 

and the governing equation of the system as, 

1 / 1 , 

mx + cx + kxlxl^^= /Qsin(cOQt + —at ) (3.26) 

In nondimensional form (with Xg^ = f^jk/va addition to the previously defined 
parameters), the governing equation becomes 

x" + 2^' + X I F 1^ = sin(ror + (3.27) 

The Simulink model for numerical solution of above equation is shovm in Appendix 5. 
Fig 3.12 shows the response of accelerating and decelerating case (r^ =0.001) and 


31 




Non-Dimensional instantaneous frequency r(= + r„r) 


Fig. 3.12 Response for nonlinear bearing system 


damping coefficient ^=0.01. The apparent shift in resonance, in this case can be seen to 
be less than that of the symmetric square case (Fig.3.8). The influence of various 
acceleration and decceleration on the response patterns are shown in Figs. 3.13 (a), (b). 
Comparison with equivalent polynomial forms (Figs. 3.14, 3.15) reveal that other 
features, including the subharmonic phenomenon also show trends similar to that of 
square symmetric nonlinearity. 
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Fig. 3.13(a) Influence of acceleration on nonlinear bearing system 



Fig. 3.13(b) Influence of deceleration on nonlinear bearing system 
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Non-Dimensional displacement x 



Fig. 3.14 CmparUon ,f bearing non-linearily | x |>^ and equivalent 

polynomial system 
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3.4 Bilinear Nonlinearity 

It is customary to represent fatigue cracks in structural members subjected to repeated 
loading through Breathing Crack models. A breathing crack is an idealisation to consider 
the change in stiffness of the structure (say, a beam), during vibration, when the crack is 
subjected to both tension and compression, during each cycle of vibration, fhe stiffness 
is less, when the crack is under tension and as a result the crack is open. The stiffness is 
higher and taken approximately equal to that of an uncracked beam, during that regime of 
the vibration cycle, when the crack is under compression and consequently closed. Such a 
system can be idealised through a bilinear governing differential equation given below. 
im^cx + kx = /oSin(<y(,t + l/2ar^) (3.29) 

A: = (1 - s)k^ fQj. X >0, 

^ = K for jc ^ 0 , 

Taking p — jm (and the previously defined parameters), and writing the equation in 
nondimensional form 

x'' + 2^' + (l-f)x =sin(r(,r + l/2r„r^), forx^O (3.30) 

x' + 2<^' + x = sin(ror + 1/2 ) for x < 0 

the response is numerically simulated through the Simulink model shown in Appendix 6. 
- Fig 3.16 shows the response for the system under acceleration for different values of 
nondimensional bilinear stiffness ((l-f) = 0.999, 0.777). With increase (1-5 ), the system 
resonance is seen to shift away from r =1.0. The figure also includes the curve under 
similar circumstances, for a purely linear case (5 =0). The presence of bilinearity (due to 
say the presence of a crack) can be seen to result in higher amplitudes at resonance, 
which is as expected due to reduction in the system stiffness. With higher values of e , 
the reduction in stiffness is more and consequently resonant amplitudes are higher. The 
response pattern, of the bilinear system, with different acceleration rates is shown in Fig. 
3.17. The response of this system is also associated with the subharmonic phenomenon 
and explanation can be attempted through equivalent polynomial forms similar to the 
previous cases. 
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Non-Dimensional displacement x Non-Dimensional displacement x 



Non-Dimensional instantaneous frequency r(= +rgT) 

Fig. 3.16 Influence of bilinear stiffness during acceleration 



Non-Dimensional instantaneous frequency r(= + r^r) 

Fig3.17 Influence of acceleration rates (bilinear system) 
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CHAPTER 4 


EXPERIMENTAL INVESTIGATIONS 

Experimental investigations have been carried out for better understanding of the 
response patterns of a nonlinear system under a sweep test. The experiments were carried 
out on a readily available thin beam, which under large amplitudes of oscillation can 
exhibit nonlinear characteristics. The exercise was carried out with the aid of Virtual 
Instruments developed in the programming software LabVlEW. 

4.1 Experimental Set-Up 

The experimental set-up comprised a thin, commercially available ruler scale mounted on 
an electrodynamic shaker at one end (Fig.4.1). The beam approximately represented a 
cantilever configuration. However, since the vibrational modes were far separated, 
analysis could be carried out for the fundamental mode by treating the response to be 
similar to that of a single degree freedom system. The geometric and material properties 
of the beam are given below. 

Length, I = 362 mm Width, w = 25 mm Thickness, d= 1.68mw 

k£ 

Modulus of Elasticity, E = 1 00 GPa Density, p = 7200 — - 


4.2 Instrumentation 

The instrumentation employed for the study is shown in Figs. 4.2. The response of the 
beam was picked up through a miniature accelerometer Type: 4374 Bruel & Kjaer, 
Charge Sensitivity: 0.153 pcfms~^. Voltage Sensitivity: 0.197 mvlms . The signal 
from the accelerometer was conditioned and amplified in a Charge Amplifier (Bruel & 
Kjaer, I'ypc 2635). I'he excitation signal for the Electrodynamic Shaker was generated in 
a Pentium PC, using UbVlEW. LabVlEW is program development application, much 
like C or BASIC. However, LabVlEW uses graphical programming language, G, to 
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Fig .4.1 Specimen mounted on electrodynamic shaker 



Fig. 4.2 Experimental set up. 
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create programs in block diagram form. In addition of having extensive libraries of 
functions, LabVIEW includes libraries for data acquisition, GPIB and serial instrument 
control, data analysis, data presentation, and data storage. The programs developed in 
LabVIEW are called virtual instruments {VI because their appearance and operation 
can imitate actual instruments, 1 hus a VJ consists of: 

Front panel. The interactive user interface of the program is called front panel, because it 
simulates the panel of a physical instrument. The front panel can contain knobs, push 
buttons, graphs, and other controls and indicators. Data can be entered using mouse or 
keyboard. 

Block Diagram: The block diagram is a pictorial solution to a programming problem. It is 
also the source code for the VI . 

The Block Diagram developed, during the course of this study, for signal generation is 
shown in Fig. 4.3 (b). The Front Panel of the VI is shown in Fig.4.3 (a). 
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Fig. 4.3(a) Front panel of VI for signal generation 



Fig. 4.3(b) Block diagram of VI for signal generation 
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The Virtual Instrument can generate a sinusoidal signal with a time varying frequency. 
Ihe initial frequency as well as the rate of change of frequency (positive or negative) can 
be set through the Front Panel of the VI. A software setting is also provided for the 
amplitude of the sinusoidal signal. The PC is equipped with an AID card (National 
Instruments, Texas, AT-MIO-E16, with a scan rate of 250000 samples/s). The signal 
generated by the VI and coming out through the CEird is taken through an accessory (BNC 
2090), which connects the A/D card to BNC type connectors, to a Power Amplifier 
(Type: SS250, MB Dynamics, Inc.) The output of the Power amplifier is fed to the 
electrodynamic shaker. The response signal is also taken to the computer, through the 
AID card. The VI for data acquisition and display is shown in Figs. 4,4 (a), (b). 
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Fig. 4.4(a) Front panel of VI used for Data Acqubition 



Fig. 4.4(b) Block diagram of VI used for Data Acquisition. 
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4.3 Experimental Results 

Prior to carrying out, nonstationary analysis, a Rap-Test was performed on the beam to 
identify its linear natural frequencies. The beam was very slightly perturbed (to keep the 
vibrations in the linear zone) at its free end and the resultant free vibrations were 
recorded. 1 he FFT of this response is shown in Fig. 4.5, from which the first two natural 
frequencies can be seen as 7.68Hz and 47.22Hz. The natural frequency were calculated 
theoretically also, employing the formula 

(4.1) 

where for a cantilever configuration =3.52, 22giving, 
a)^ =7.72Hz 
Gij =48.29Hz 

The theoretical and experimental results can be seen to be very close. Also the natural 
frequencies can be seen to be far separated. 

In order to confirm that the beam behaves as a nonlinear system for large amplitude 
oscillations, it was excited by a harmonic force of frequency 10 Hz and relatively large 
amplitude as shown in Fig.4.6 (a). The response FFT shown in Fig.4.6 (b) shows a 3m 
harmonic at 30 Hz., in addition to the harmonics at 20 Hz and 40 Hz., corresponding to 
the excitation frequency. Such test was carried out for various excitation frequencies, and 
the FFT of the response was characterised by the presence of 3m harmonic, confirming 
the nonlinear behaviour of beam at large amplitude vibrations. 

The sweep tests were then carried out for different acceleration and deceleration rates. 
Fig. 4.7(a) shows the excitation function, generated in the signal generator. The excitation 
frequency is swept from 0 Hz to 10 Hz. at a constant rate in 60 seconds. The response is 
shown in Fig. 4.7 (b). The peak vibrations can be seen to occur before the linear natural 
frequency (7.68 Hz). I'he excitation for deceleration from 10 Hz to 0 Hz in 60 seconds, is 
shown in Fig. 4.8 (a) and the corresponding response is given in Fig. 4.8 (b). The peak 
response while, decelerating can be observed to be at an instantaneous frequency lower. 
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Fig. 4.5 Rap test results for identiflcation of natural frequencies of beam. 


than that while acceleration. Also, while deceleration the peak amplitude is higher than 
that during acceleration. This indicates a behaviour similar to a soft spring, with 
backbone curve bent to the left. The jump phenomenon can also be noticed from the 
figures. Figs. 4.7 (b) and 4.8 (b) also show the 1/3 subharmonics very clearly. The 14 
subharmonic, if present seems to be very weak. This suggests that the system can be 
approximately modeled through an equation of the type 
nvi + cx + k^x-kyX^ = J\ sinaX 

While -yJk^Tm = <y, , the first linear natural frequency of the system has already been 
identified during the Rap- l est as 7.68 Hz, no attempt has been made in the present study 
to identify the nonlinear stiffness value k ^ . 
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Fig. 4.6(b) Response exhibiting nonlinearity of beam 
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Fig. 4.7(b) Noil-stationary response of beam 
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Instantaneous frequency Hz 


Fig. 4.8(a) Excitation force with frequency swept from 10 to 0 Hz in 60 seconds 



0123456789 10 

Instantaneous frequency Hz 


Fig. 4.8(b) Non-stationary response of beam 
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acement 




Superposisition of the response during acceleration and deceleration, as shown in Fig. 
4.9, may help in providing information on the backbone curve. However, construction of 
such a curve from the experimental data, has also not been attempted during this study. 

The response was obtained for different acceleration and deceleration rates. Fig. 4.10 (a), 
(b), show the response for acceleration and deceleration respectively, between 0 and 10 
Hz in 100 seconds. Figs. 4.1 1 (a), (b) show the corresponding behaviour, when the transit 
time was reduced to 60 seconds while Figs. 4.12 (a), (b) show the response behaviour 
when transit time was further reduced to 40 seconds. It can be seen, through a 
comparative analysis of Figs. 4.8, 4.10 to 4.12, that with faster transit rates the peak 
response lowers in magnitude ih addition to showing greater shift away from the linear 
natural frequency value. 
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Instantaneous Displacement Instantaneous Displacement 
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- 0.8 
-1 

Instantaneous frequency Hz 

Fig. 4.10(a) Non-stationary response during acceleration 
(0 - 10 Hz in 100 seconds) 



Instantaneous frequency Hz 

Fig. 4.10(b) Non-stationary response during deceleration 
(10 - 0 Hz in 100 seconds) 
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Instantaneous Displacement Instantaneous Displacement 



Instantaneous frequency Hz 

Fig. 4.11(a) Non-stationary response during acceleration 
(0 -10 Hz in 60 seconds) 



Instantaneous frequency Hz 

Fig. 4.11(b) Non-stationary response during deceleration 
(10 — 0 Hz in 60 seconds) 
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Instantaneous Displacement Instantaneous Displacemei 



6 6.5 7 7.5 8 8.5 9 9.5 10 

Instantaneous frequency Hz 

Fig. 4.12(a) Non-stationary response during acceleration 
(0 -10 Hz in 40 seconds) 



Instantaneous frequency Hz 

Fig. 4.12(b) Non-stationary response during deceleration 
(10 -0 Hz in 40 seconds) 
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CHAPTER 5 


CONCLUSIONS 

The objective of the present study was to carry out numerical and experimental exercises 
to understand the behaviour of nonlinear systems subjected to harmonic excitation with 
constant variation in frequency. Some commonly encountered nonlinearities in 
mechanical systems were considered. An attempt was made to understand their behaviour 
by comparison with that of linear systems. The polynomial form of nonlinearity, 
specifically the cubic nonlinearity of a Duffing oscillator was further treated as 
benchmark for understanding peculiarities of other nonlinear systems. It was observed 
that the nonstationary behaviour of nonlinear systems is different from linear systems, in 
a variety of ways. A major feature of the nonlinear systems is that the response patterns, 
during acceleration and deceleration are not symmetric. The degree of asymmetry, is 
dependent on the type of nonlinearity, in addition to the rate of frequency change. For 
polynomial forms of nonlinearity, the assymmetry of these patterns is governed by the 
backbone curve. Lack of analytical models for other forms of nonlinearities, restricts their 
analysis to attempting piecewise polynomial approximations. These piecewise 
polynomial models help in justifying the presence of subharmonics in the response. 

Experimental investigations, in addition to nonstationary vibration analysis, were also 
aimed at providing exercise in Virtual Instrumentation for Signal Generation as well as 
Data Acquisition. This task was successfully accomplished. The quality of the 
experimental results was observed to be good as distinct features of nonlinear behaviour 
during sweep tests could be clearly observed and rationalised. 

Future numerical exercises can be carried out by incorporating damping nonlinearities 
also, in the models. The simplicity of the MATLAB Simulink models, make this task 
easy. Experimental investigations, in further studies can focus on systems like bearings, 
shafts etc. The Virtual Instruments developed can also be employed to investigate 
influence of cracks of various sizes and locations, on the nonlinear dynamics of beams. 
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APPENDIX 1 



SIMULINK MODEL FOR LINEAR SYSTEM 





APPENDIX 2 



SIMULINK MODEL FOR CUBIC NONLINE. 







APPENDIX 3 



SIMULINK MODEL FOR SQUARE-CUBIC NONLINE 








APPENDIX 4 



SIMULINK MODEL FOR SQUARE SYMMETRIC NONLINEARITY 





APPENDIX 5 



SIMULINK MODEL FOR BEARING NON L. 







APPENDIX 6 



SIMULINK MODEL FOR BILINEAR NONLINEARITY 





